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We consider the entanglement properties of the quantum phase transition in the single-mode 
superradiance model, involving the interaction of a boson mode and an ensemble of atoms. For 
infinite system size, the atom-field entanglement of formation diverges logarithmically with the 
correlation length exponent. Using a continuous variable representation, we compare this to the 
divergence of the entropy in conformal field theories, and derive an exact expression for the scaled 
concurrence and the cusp-like non-analyticity of the momentum squeezing. 
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Entanglement has shot to prominence in recent years 
on the back of the success of three key areas: quantum 
computing, quantum cryptography, and quantum tele- 
portation. In this quantum information paradigm, entan- 
glement is a resource which can be exploited to perform 
hitherto unimagined physical tasks. 

Latterly, a new emphasis has emerged in which en- 
tanglement is related to properties of interacting many- 
body systems. This approach is being pursued most 
vigourously in connection with quantum phase transi- 
tions (QPTs) Q, as it is hoped that entanglement may 
shed light upon the dramatic effects occurring in criti- 
cal systems which, by their very nature, involve complex 
collective quantum mechanical behaviour. A complete 
theory of many-body entanglement is still lacking. Cur- 
rent techniques are reliant upon bipartite decompositions 
of the total system, and the criteria for selecting the most 
pertinent decomposition are by no means clear. 

Investigations so far have therefore been restricted to 
interacting spin-1/2 systems on a one-dimensional lattice 
0i Hi S U or on a simplex |^ , which require the (more 
or less artificial) splitting into two spin-subsystems. 

In this Letter, we study the entanglementproperties of 
the one-mode superradiance (Dicke) model |3 , where col- 
lective and coherent behaviour of pseudo-spins (atoms) 
is induced by coupling (with interaction constant A) to 
a physically distinct single-boson subsystem. We present 
here exact solutions for the entanglement of formation 
between these two subsystems, and for the pairwise en- 
tanglement between atoms at and away from the criti- 
cal point Ac. Recently the QPT in this model has been 
related to the emergence of chaos for A > Ac in a corre- 
sponding classical Hamiltonian Our real-space rep- 
resentation of the modes allows us to analyse the scaling 
of the atom-field entanglement at the critical point, and 
to compare with results from conformal field theories for 
one-dimensional spin chains 0. Furthermore, we derive 
explicit expressions for the concurrence and the related 
(momentum) squeezing for all coupling parameters A. 

A model that has drawn considerable interest in the 
context of entanglement near criticality is the XY-mode\. 
In ferromagnetic spin 1/2 chains, the concurrence as a 
function of system size been used [3| to demonstrate scal- 



ing of entanglement near the transition point. Osterloh 
et al. |3] have shown that the derivatives of the con- 
currence between neighbour and next-nearest neighbour 
spins exhibits a universal scaling behaviour in the re- 
gion of the critical point in this model. Furthermore, the 
study of such systems has led Osbourne and Nielsen Q 
to the notion of a 'critically entangled' system where the 
correlation length ^ of the system is divergent and entan- 
glement exists over all length scales. Vidal et al. |i| have 
used an alternative approach and studied the entangle- 
ment between blocks of L contiguous spins and the rest 
of the chain and have found a striking relation to the en- 
tropy Sl ~ (c + c) /6 log L + const in I -I- 1 conformal field 
theories with central charges c and c. 

We start by describing our model, which is the single- 
mode Dicke Hamiltonian describing the interaction of N 
two-level atoms of splitting with a single bosonic mode 
of frequency uj 
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where the second form follows from the introduction of 
collective spin operators of length j — N/2. In the ther- 
modynamic limit, N,i oo, the system undergoes a 
QPT at a critical coupling of A = Ac = ^/aJu^/2, at which 
point the system changes from a largely unexcited nor- 
mal phase to a super-radiant one in which both the field 
and atomic collection acquire macroscopic occupations. 

Similar to the large-spin problem analysed in this con- 
text Q , the Dicke Hamiltonian can be regarded as a zero- 
dimensional field theory with mean-field type behaviour, 
where the Sn permutation symmetry of the atoms and 
the absence of an intrinsic length scale makes the model 
exactly solvable. Despite this simplicity, our model ex- 
hibits many non-trivial properties; in particular, exact 
solutions for the non-analyticities of the entanglement of 
formation and the concurrence can be related to the scal- 
ing exponent, the finite-size behavior, and the underlying 
semi-classical integrable/chaos cross-over which has been 
shown to occur around the phase transition 0. 
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FIG. 1: Entanglement of formation Soo between atoms and 
field for both N oo and finite A''. Inset: Scaling of the 
value of the entanglement maximum as a function of logj A'^. 
The Hamiltonian is on scaled resonance lj — too = 1. 

The starting point for our analysis in the thermo- 
dynamic limit is the HoIstein-PrimakofF representation 
|9| of the angular momentum operators Jz = (6^6 — j), 
J+ = \/2j — b^b, J_ = j|. Here, b and b^ are bosonic 
operators that convert Ti. into a two-mode bosonic prob- 
lem. This allows us to obtain effective Hamiltonians 
that are exact in the thermodynamic limit, by neglecting 
terms from expansions of the Holstein-Primakoff square- 
rootS/S-l. In the normal phase, A < Ac, we expand the 
square-roots directly and obtain the effective Hamilto- 
nian 

7^(1) = ujQb^b + uja^a + A (a^ + a) (6^ + 6) - jujo. (2) 

In the super-radiant phase, we first displace both boson 
modes by quantities proportional to ^/j before we ap- 
proximate the square-roots. This leads to a second effec- 
tive Hamiltonian, the form of which is also bilinear and 
similar to Eq. 

We now consider the normal phase ground state in 
some detail; the super-radiant phase results following 
with slight modification. The eigenstates of Ti'^^^ are 
two-mode squeezed states. Via the introduction of a 
position-momentum representation for the two oscilla- 
tors; X = {a^ + q) , y = ^2ui (^^ ^) ' ^it^ 
menta defined canonically, we may write the ground-state 
wavefunction as 

where = ^ + ujq ± y/ (wq — lo^)^ + IGX'^ujuJoj are 
the excitation energies of the system, s = sin 7, c = cos 7, 

and the angle 27 = arctan[4A^a;a;o/(wo ~ '^^)] charac- 
terises the squeezing axis. This wavefunction forms the 
basis of the current analysis. 



Entanglement of formation - As a measure of the en- 
tanglement between the atoms and the field, we calculate 
the von-Neumann entropy S = —tip logj p of the reduced 
density matrix (RDM) p of the field-mode. In the nor- 
mal phase, p is simply determined by the ground state 
wave function, Eq. (|3J), whereas in the super-radiant 
phase two degenerate ground states exist that have wave 
functions ^± similar to Eq. Q, but displaced from 
the origin by amounts proportional to ±-\/7. This de- 
generacy arises from the breaking of the parity symme- 
try H = exp {in [a^ a + Jz + j]} for A > Ac- Because 
'i'^ and '^^ are orthogonal, the convex nature of the 
von-Neumann entropy [l(j| implies that in the SR phase 
S{pcat) = S{p±) + Ij where p± is the RDM of either 
of the two (macroscopically separated for large N) so- 
lutions, and Pent is the RDM of the superposition 'cat' 
state of the two. The cat state restores the broken parity, 
and thus the latter expression will be used for comparison 
with the numerical results for finite N. 

Having clarified this additional distinction between the 
two phases, we now explicitly calculate the normal phase 
RDM in the x-representation, 

/oo 
dyfLivW {x,y)^{x',y). (4) 
-00 

Here, cl is a normalisation constant, and the introduc- 
tion of the cut-off function = will allow us 
to discuss the effect of a partial trace over the atomic (y) 
modes (see below). A straightforward calculation shows 
that PL is identical to the density matrix of a single har- 
monic oscillator with frequency 17^ in a canonical ensem- 
ble at temperature T = where 

cosh/317L = 1 + 2 ,y ,\ ^' — . 5 

(£_ — £-|_)'^C^S'^ 

The entropy obtained from is thus given by the 
expression 

SdO = [CcothC-ln(2sinhC)]/ln2, C = P^l/2. (6) 

This strikingly simple result allows some interesting ob- 
servations. First of all, the entropy Soo (cut-off L — 00) 
undergoes a divergence at the critical point as we ap- 
proach Ac from either side. In the region near Ac, 
the excitation energy £_ vanishes as £_ cx |A — Acp", 
with the exponent = 1/4 describing the divergence 

— 1/2 

of the characteristic length ^ = £_ . Using 6*00 (C) = 
[I-ln(2C)+CV6]/ln2+0(C4) and C = y/^[l+0{eoo)] 
with £00 = 2£_/(£+s^c^), we find that Soo diverges loga- 
rithmically as 5*00 oc — (1/2) log2(2£oo) and hence (omit- 
ting constants), 

S'oo oc -i/log2 |A - Ac| = logaC, = 1/4. (7) 

Thus, the entanglement between the atoms and field di- 
verges with the same critical exponent as the character- 
istic length - a clear demonstration of critical entangle- 
ment. 
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As we approach Ac the parameter C = hVLoc/kBT of 
the fictitious thermal oscillator approaches zero, indicat- 
ing that a classical limit of the field RDM is being ap- 
proached, interpreted either as the temperature T go- 
ing to infinity, or the frequency Q.00 approaching zero. 
In terms of the original parameters of the system, the 
dependence of the entropy is through the ratio of ener- 
gies £00 oc e_/e+0|. Although the entanglement calcu- 
lated here is a genuine quantum property of the combined 
atom-field system, the above results highlight that in the 
limit of A'^ ^ 00 atoms, the exact mapping of the sys- 
tem to two coupled oscillators leads to emergent pseudo- 
classical behaviour. This is consistent with the observa- 
tion that the Dicke Hamiltonian in fact is strongly con- 
nected to a classical (cusp) singularity in the catastophe 
theory sense [l^ . 

We next compare the analytical result from Ea. (|5l6() 
for the entropy Soo (corresponding to completely tracing 
out the atomic mode) with the corresponding finite A^ 
results obtained from numerical diagonalisation. Fig. 1 
shows these results, and illustrates the finite size scaling. 
Defining A'^ as the position of the entropy maximum, 
and as the value of the maximum entropy, we observe 
A'^ - Ac cx iV-"-75±o.i^ and Sm oc (0.14 ± 0.01) log2 A^. 

The accuracy of the exponents are limited by the avail- 
able numerical data. The divergence of the entropy is 
logarithmic due to the symmetric nature of the spin sys- 
tem. The entropy here saturates at a maximum value of 
log2(A'^ + l)j in contrast with general spin sytems which 
saturate at log2(2^) due to their larger Hilbert spaces. 
This distinction is expected to be important in determin- 
ing the complexity of classically simulating a quantum 
phase transition fl^ 0| . An explicit plot of the en- 
tropy scaling is shown as an inset in Fig. 1, while the 
scaling of the position of the maximum point is shown as 
an inset in Fig. 2. 

We next consider the system at the critical point but 
keep the 'tracing parameter' L finite. This corresponds to 
a situation where the trace over the (atomic) y-coordinate 
is performed over only a finite gaussian effective region of 
size L for the atomic wave function. With e_ = 0, the rel- 
evant dimensionless energy scale is now £l = 2/(i^£+c^), 
and the entanglement entropy diverges as (again omitting 
constants) 

Sl (X -(1/2) log2(2ei) = log2 L, L ^ 00. (8) 

This result can now be compared with a recent calcu- 
lation by Vidal et al. of the critical entanglement 
of formation of blocks of L spins in one-dimensional in- 
teracting XY and XXZ spin-chain models. There, the 
prefactor for the logi dependence of Sl at criticahty 
is given by the central charges of the underlying con- 
formal field theory in 1 -I- 1 dimensions. Note, however, 
that a direct comparison would require the tracing out 
of L atoms from the A^-atom Hamiltonian (see below) 
with N —f 00, L fixed, but the general principle is the 
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FIG. 2: Scaled pairwise concurrence Cn = NC between two 
spins for both N ^ oo and finite N. Inset: Scaling of the 
value {+) and position (x) of the concurrence maximum, and 
the position of the entropy maximum (*) as functions of N. 
The Hamiltonian is on scaled resonance ui = u>o ~ 1- 

same. In this context, the Dicke model corresponds to a 
zero-dimensional field theory and, for N oo, is in fact 
closely related to Srednicki's simple two-oscillator rnodel 
in his introductory discussion of entropy and area [llj . 

Pairwise entanglement and concurrence.- To observe 
the behaviour of the entanglement between atoms within 
the ensemble, we proceed by considering the 'pairwise' 
entanglement of formation for mixed states, as parame- 
terized by the concurrence |l5j . The absence of an intrin- 
sic length scale in our model simplifies our calculations, 
enabling us to employ the prescription set out for sym- 
metric Dicke states in jl6| . The matrix elements of the re- 
duced density matrix pi2 for any two atoms is determined 
by the expectation values of the collective operators, ( Jz) , 
(J^), and (J^). We then define the scaled concurrence 
as Cn = NC, with C = max{0, Ai — A2 — A3 — A4}, 
where the A.^ are the square roots of the eigenvalues (in 
descending order) of pi2{criy ® o'2y)Pi2(o'iy ® ^^y)- Wang 
and Sanders have shown that Cn can be expressed in 
terms of the Kitagawa-Ueda spin squeezing |17| for sym- 
metric multi-spin states 18]. 

We show numerical results for the scaled concurrence 
Cn in Fig. 2, together with the analytic thermodynamic 
limit result described below. For all A and A^, Cn is 
less than that of the pure H^-state |j, ±(j — 1)), which 
has Cn = 2^ the maximum pairwise concurrence of any 
Dicke state [ig. For small coupling A, we recognise an 
A^-independent behaviour of Cn which may be obtained 
from perturbation theory in A as 

CN(A->0)-2aV(l + «^), a = \/{ijj + ojQ). (9) 

As with the entropy, we can perform a finite scaling anal- 
ysis of the numerical data. Again, two power law expres- 
sions are found for A'^ and C^; A"^ - Ac oc iV-o.68±o.i 
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andCM(Ac)- 



Cm oc Ar-0-25±o oi. Plots of this behaviour 



are shown as an inset in Fig. 2. 

In the thermodynamic limit N 
currence can be expressed as 



oo, the scaled con- 



C' 



{i + f,)[({d^f)~{dU)] + -ii-f,), (10) 



setting ^ = 1 and = in the normal phase (A < Ac), 
and = (Ac/A)2 and rft ^ fet + ^Ar(l - ^)/2 in the SR 
phase (A > Ac). Recalling b^ = ^wo/2(j/ — ipy/ujo), we 
can further transform Eg. 1)10(1 to establish a relation be- 
tween the scaled concurrence, the momentum squeezing 
{Apyf = (pD - {py)\ We have 



to quantify entanglement across a quantum phase tran- 
sition. The clear physical distinction between the sub- 
systems (pseudo-spin or two-level system and bosonic 
mode) enables us to see distinctly the logarithmic di- 
vergence of the entropy in the thermodynamic limit as 
a function of the coupling constant. We mention that 
quantum phase transitions have also been discussed very 
recently in the context of entanglement generation (e.g., 
for atoms in optical lattices [20]), and quantum compu- 
tation schemes A further though still mainly unex- 
plored aspect remains the fundamental role of the phase 
transition for the connection between entanglement and 
the underlying integrable to quantum chaotic transition. 



Co 



(1 + m) 



{Apyf/UO 



+ -(1-m), (11) 



where again, setting ^ = (Ac/A)^ gives the super-radiant 
phase equivalent. The concurrence can be explicitely re- 
lated to the parameters in the reduced y (atom) oscil- 
lator density matrix as Coo = 1 ^ (/iil/tJo) coth(/?il/2) 
with cosh/317 = 1 + 2e-£+/D, D = [cs(e_ - e+)f, 
and 20,/ sinh/30 = D/{e-C^ + s+s^)- Due to symme- 
try, these are the same parameters as for the reduced 
field {x) density matrix poo, Eq. with s = sin 7 and 
c = cos 7 interchanged. After simple algebra one obtains 
Coo = 1 — M(e-s^ + e+(?)/uJo. Due to space restrictions, 
we only give analytical results at resonance {uj — loq)^ 



CX<1 
00 



X = A/Ac (12) 



V2x2 



(sin2 7) \j I + x^ ~ \f{l - x^f 



(cos^ 7) \jl + x'^ + \J{l-x^f 



(13) 



where 27 = arctan[2/(a;2 — 1)] in the SR phase. These ex- 
plicit expressions reveal the square-root non-analyticity 
of the scaled concurrence near the critical point Ac- The 
concurrence assumes its maximum Coo = 1 ~ V2/2 ~ 
0.293 at the critical point A = Ac. We note that 
Eq. ((12(1 is consistent with the maximum of the (un- 
sealed) concurrence approaching the critical point in a 
related, dissipative version of the Dicke model in the nor- 
mal Dhase|l9j. Our findings are also in agreement with 
the behaviour of the concurrence in the collective spin- 
model, H = ^{2\/N){Sl+-iSl)-2S^ + {\/2){l + -i) |i| 
and differ from ID spin chains, where the maximum of 
the C does not coincide with its non-analyticity at the 
critical point. We also note here that the squeezing ob- 
tains its minimal value at Ac, which is again in agreement 
with the above spin model. 

In conclusion, we have obtained exact results for the 
entropy and the concurrence in a model that allows us 
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